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1. Introduction 

It is well known that the transition from a primordial inflationary phase to a 
decelerated one, typical of our present cosmological evolution, is associated with 
the production of a cosmic background of relic gravity waves-^"^. The spectral 
distribution of their energy density may provide direct information on the very 
early history of our Universe, and can be used, in particular, to reconstruct the 
time dependence of the Hubble parameter^. 

Deflation, however, is not the only violent process typical of primordial evo- 
lution able to amplify a metric fluctuation. Although less known (or less studied, 
up to now at least), it is a fact that gravitons can be produced from the vacuum 
also as a consequence of a phase of dynamical dimensional reduction®'^, in which 
a given number of "internal" dimensions shrink down to a final compactification 
scale. Another possible process which may lead to a cosmological graviton produc- 
tion, and which (to our knowledge) has not yet been discussed in the literature, is 
the time variation of the effective gravitational coupling constant G. 

The main purpose of this paper is to compute the expected spectrum of the 
cosmic gravitons background, by including both the contributions of dimensional 
reduction and of G among the possible sources (besides inflation), and by using 
a Brans-Dicke-like graviton-dilaton coupling as a dynamical model of variable G. 
We are led to this choice, in particular, by the models of early universe evolu- 
tion based on the low energy string effective action which suggest that the 
standard radiation-dominated cosmology is preceeded by a dual, "string- driven" 
phase, in which the effective gravitational coupling changes just because of the 
time-dependence of the dilaton background. The possibility of looking for tracks 
of such a string phase in the properties of the cosmic graviton spectrum provides 
indeed one of the main motivations of the present work. 

The paper is organized as follows. In Section II we deduce the linearized 
equation for a gravitational wave perturbation in a Brans-Dicke background, and, 
in Section III, this equation is used to compute the spectral distribution of the 
gravitons, produced by the cosmological background transitions. We shall take 
into account the dilaton-driven variation of G in a higher dimensional framework 
in which also the scale of the internal spatial dimensions is allowed to vary, and 
in which the matter-dominated and radiation-dominated evolution of the external 
space follows a phase of accelerated (i.e. inflationary) expansion. The squeezing 
parameter^^ corresponding to this scenario will be given in Section IV. 

In Section V the present bounds on the energy density distribution of the relic 
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gravitons are used to obtain information, and constraints, on the value of the curva- 
ture scale at the transition between the inflationary and the radiation-dominated 
era, versus the parameters characterizing the background kinematics. The pre- 
dictions of some string-inspired cosmological models (and of related Kaluza-Klein 
scenarios) will be compared with these bounds in Section VI. The main conclusions 
of this paper will be finally summarized in Section VII. 



2. Gravitational perturbations of a Brans-Dicke background. 

The starting point to discuss the production of gravitons, induced by a cosmo- 
logical background transition, is the linearized wave equation for a gravitational 
perturbation propagating freely in the given background. In order to include the 
effect of a changing gravitational coupling, such an equation will be obtained by 
perturbing (at fixed sources) the Brans-Dicke field equations, around a background 
configuration which includes a time-dependent dilaton field. 

It should be perhaps recalled that, in a general relativity context and in 
a spatially fiat Priedmann-Robertson- Walker manifold, a gravity-wave perturba- 
tion obeys the same equation as a minimally coupled massless scalar field^'^^'^^. 
In a Brans-Dicke context, however, the graviton wave equation is different from 
the covariant Klein-Gordon equation, because the gravitational perturbations are 
coupled not only to the background metric tensor, but also to the scalar dilaton 
background (/)(t) representing the G variation. 

Our background field dynamics is assumed to be described, in D dimensions, 
by the following scalar-tensor action 

S = J d^'x e-'^ ^MiR - ^g'^'d^.d,^) + Sm (2.1) 

where u) is the usual Brans-Dicke parameter, and Sm represent the possible con- 
tribution of matter sources, with ^/\g\TfJ_^ = 25Sm/Sg^'^. The variation of this 
action with respect to provide the dilaton equation 

R + io{V(l)f -2ion(j) = (2.2) 

where V denotes the Riemann covariant derivative, and □ = g^'^V The vari- 
ation with respect to g^i,, combined with (2.2), provides the equation 

i?M " + V^V^c^ + (a; + 1)K((V0)2 - □</.) - V^0V>] = MGne^T^ (2.3) 
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Note that we are using an exponential parametrization for the dilaton field, to 
make contact with the string cosmology models. For a; = oo, = const we recover 
general relativity, while for a; = — 1 eq (2.1) reduces indeed to the truncated low 
energy string efi'ective action with phenomenological matter sources^°~^^. 

The free linearized wave equation for a metric fluctuation Sg/^i, = h^i, is now 
obtained by perturbing eqs.(2.2) and (2.3), keeping all sources (dilaton included) 
fixed, 5T^ = = 5(j). It should be stressed that we have not explicitly included in 
the action a possible dilaton potential term, F(0), as its contribution to the per- 
turbation is vanishing for 50 = 0. It is true that, in a class of duality-symmetric 
string cosmological models^°~^^'^^, the dilaton self interactions may also occur 
trough a coupling to the metric, and lead to a two loop potential of the form 
V = Vo[exp (20 — 21n y^)], for which 5V oc V5g. This potential, however, is 
expected to affect in a significant way only the transition region between the infla- 
tionary and the radiation dominated regime^^. Its contribution to the perturbation 
equations may then be neglected for the purpose of this paper where, as discussed 
in the following section, we shall evaluate the graviton spectrum in the "sudden" 
approximation, namely in the approximation in which the kinematic details of 
the transition regime are ignored, and the rapid exponential decay of the high 
frequency tail of the spectrum is replaced by a suitable high frequency cutoff. 

We perform then the transformation g^^i, — > g^,^ + 5gfj,i,, with 

5g^,u^K. , d4>^0^ST^'' (2.4) 

By neglecting corrections of order higher than first in /i^j^ (so that, for instance, 
Sg^'^ = —h'^^), we are led to the following variational expressions 

5{ucj>) = -h^-v^v,(j> - ^'^"a«#r^, 

5{V^V^4>) = -/i"" V^Va0 - g^'^dp<j>5Tl^ 

(5(V^0V» = -h^'^d^^da^ (2.5) 

Here 

5r^. = Ig'^^i^^Kp + V. V - V/3 V) (2.6) 

and 5R^i, is the linearized expression {oiR^jy{5g) (note that all covariant deriva- 
tives, as well as all operations of raising index on h^^, are now to be understood 
as performed with the help of the background metric g,j,y). 
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We choose, in particular, a time dependent background with cj) = (f){t), and 
a homogeneous diagonal metric describing a general situation of dimensional de- 
coupling , in which d dimensions expand with scale factor a{t), and n dimensions 
contract with scale factor b{t). In a synchronous frame 

900 = '^ , 9ij = -a'^{thij{x) , 9ab = -b'^{thab{y) 

9o^.^^^9ia , (P = (p{t) (2.7) 

(conventions: = 1, D = d+n + 1; i, j = 1, d; a,b = 1, n; t is the cosmic 
time coordinate, and 7^, ^ab are the metric tensors of two maximally symmetric 
euclidean manifolds, parametrized respectively by "internal" and "external" coor- 
dinates {x*} and {2/"}). We are interested moreover, in a pure tensor gravitational 
perturbation, decoupled from sources, representing a gravitational wave propagat- 
ing in the d-dimensional external space, such that /i^i/ = h^y{x, t), Hq^ = = hafj,, 
and which satisfies the transverse, traceless gauge condition 

5^"V = = V./i^" (2.8) 
In this case we have, for the background (2.7), 

(5i? = , 5V% = ~h^, (2.9) 

(a dot denotes derivative with respect to t). The perturbation of eq.(2.2) is thus 
trivially satisfied, while the perturbation of eq.(2.3) provides for /i^j, the linearized 
wave equation 

5R^ + ^(/>V^7"° - /i"" V^V«(/> = (2.10) 

which, being u independent, is remarkably the same for all Brans- Dicke models. 

The non- vanishing components of the background Ricci tensor, for the metric 
(2.7), are given by 

Rq ° = -d{H + H'^) - n{F + F^) 
Ri^ = -^Ri^ {-fix)) - 5lidH'^ + H + nHF) 

Ra ' = -^Ra Hl{y)) - S'ainF^ + F + dHF) (2.11) 

where H = a/a, F = h/h and -R(7) denotes the Ricci tensor for the n dimensional 
euclidean spaces computed from the metrics jij{x) and ^ab{y)- By using the 
relations 

gij=2Hg,j , g'^ = -2Hg'^ (2.12) 
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one obtains-'^^'-'^^, to the first order in Sgij = hij, 

5{g^) = Sig^'^gik) = = {g^^hik)- (2.13) 
It is thus simple to show (in the gauge g^^hij = ) that 

^i? = ^ 5{g'^9ik) = = 5H 
SH' = ^ Sig^'gikf = 

S{H''Si^) = ^Hhi^ (2.14) 

(the corresponding perturbations of the F terms are aU vanishing, since Sgab = 0). 
Therefore 

= = SRa ^ 

SR,^ = -d{^) - ^Hhl - - "^Fk^ (2.15) 

We shall consider, in particular, a flat euclidean metric 7jfe = Sik, so that 
T^j{x) = = Ri^{'^). The gauge condition W{'j)hi^ = reduces to djhi^ = 0, 
and implies^^'^^: 

SRi^ = -^V^hi^ (2.16) 
with = S^Widj. We thus recover the usual result 



6R,^ = -l[hi^ + {dH + nF)hi^ - \v^h,'] = -^n/ii^' (2.17) 



valid whenever the background is isotropic in the polarization plane, orthogonal 
to the direction of propagation of the wave-*^*. 

On the other hand we have, for the background (2.7), 

ViV,-0 = ]^^gij (2.18) 

Moreover, by using eq.(2.12), 

g^^hik-y^gih^hi^ (2.19) 
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The linearised wave equation (2.10) thus reduces to 



nhi^ -^hi^ = (2.20) 
and, in terms of the eigenstates of the Laplace operator, 

V^hi^{k) ^ -k^hi^{k) (2.21) 

it takes the form 

hi ^ + (dH + nF- ^)hi ^ + i-fhi ^ = (2.22) 

a 

For later applications, it is convenient to rewrite this equation in term of the 
conformal time coordinate r/, defined by dt/dr] = a. Denoting with a prime the 
differentiation with respect to rj, and defining 

V;,^=/ti^a('^-i)/26"/2e-'^/2 (2.23) 

we get finally, from eq.(2.22), that each polarization mode i/jj (k) must satisfy the 
equation: 

^" + (/c^ - V)ij = (2.24) 

where 

-„(„-2)(-)^ + -^' +-„(d-l)---(d-l)-^'---^' (2.25) 

This effective potential generalizes to a higher number of dimensions the four- 
dimensional equation, used by Grishchuk and coUaborators^'"^'^, to study the cos- 
mological amplification of the quantum fiuctuations of the metric tensor. In addi- 
tion, it takes into account the coupling of the metric perturbations to a possible 
time variation of the gravitational coupling constant (</>' 7^ ), and to a possi- 
ble variation of the scale of n "internal" compactified dimensions {b' 7^ 0). It 
may be interesting to note that this potential can also be expressed in terms of 
the scale factors only, by eliminating the explicit dilaton dependence through the 
background equation (2.2), which implies 

2 4^ 2^ ^ 2 b^ 



— [2d— + 2n^ + did - 3)(-)2 + n{n - iM^^f + 2nid - 1)^1 (2.26) 
A.ijj a a ah 

In this way one can re-introduce the w-dependcnce which is otherwise hidden in 

the particular choice of the dilaton background. For the purpose of this paper, 

however, it will be more convenient to work directly with the form (2.25) of the 

potential, in which cj) appears explicitly. 



3. Parametrization of the graviton spectrum 
for a general model of background evolution 

As discussed in the previous section, the present day background of cosmic 
gravitational waves may include, among its sources, not only a metric transition 
(deflation, dynamical dimensional reduction), but also a dilaton transition between 
two (or more) regimes with different gravitational coupling. 

In order to take all these contribution into account, we shall consider the 
background metric of eq.(2.7) (with flat maximally symmetric subspaces 7^^- = 
^iji lab = Sab), Starting with an initial configuration in which, for t] < — d 
dimensions inflate with scale factor a{r]), n dimensions shrink with scale factor 
b{r]), and the dilaton coupling is growing according to 

a r]~'^ , h ^ r]^ , (p ^Ina , r/ < — ryi (3-1) 

(note that in this equation rj ranges over negative values, so that a, (3 and 7 
are all positive). We shall assume that this phase is followed, at 77 = —771 and 
rj = rj2 respectively, by the standard radiation-dominated and matter-dominated 
expansion of three spatial dimensions. During these two last epochs, however, the 
gravitational coupling and the compactification scale of the possible additional ni 
internal dimensions are not assumed to be frozen, but they are allowed to vary as 

a ~ 7/ , h ^ r]~^^ , ~ 71 In a , — 771 < 77 < 772 

a rf' , 6 ~ 'q~^'^ , ~ 72 In a , < 772 < 77 (3-2) 

According to this model of background evolution, the effective potential (2.25) 
becomes 

^(^) = [["(^ - 1 - 7) - n/3 + 1]2 - 1] , r/ < -771 
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^(?7) = + 71 -1)^-1]] , -Vi<V<V2 

V{v) = ^[{niPi + 2^2-Sf-l] , V2<V (3.3) 

(note that it goes to zero as 77 ^ ±00). A particular solution of eq.(2.24) for ^p{k) 
can thus be written in terms of the first and the second kind Hankel functions 
H^^^ and if*^^^ (we follow the notation of Ref.(19)), ifj^k^rj) ~ r}^ H^'^\krf), which 
correspond to free oscillating modes in the |?7| 00 limit, as t]^ H^'^'^\kri) — > 
e^^^"^ /y/k (the minus and plus sign corresponds, respectively, to H^"^^ and H^^^). 

The effective potential barrier (3.3) leads to an amplification of the gravita- 
tional perturbations or, equivalently, to a graviton production from the 
vacuum^'^'^"^. Indeed, starting with incoming modes which are of positive fre- 
quency with respect to the vacuum at the left of the barrier (rj — > —00), one has 
in general, for rj — > -|-oo, a linear combination of modes which are of positive and 
negative frequency, with respect to the vacuum at the right of the barrier. The 
superposition coefficients c±{k) define the Bogoliubov transformation^^ connect- 
ing the "left" and "right" vacuum, and determine the spectral distribution of the 
produced gravitons. 

By assuming, in our case, the " in" states of the gravitational field correspond 
to the Bunch-Davics "conformal" vacuum^'^'^'^, we can write the general solution 
of eq.(2.24), for each mode ipi^): the three temporal regions as follows: 

'il;i{k) = Crr^Hl^^ , v<-m 

iPnik) = ii^A+H^^\krj) + A_Hl'\krj)] , -m < V < V2 

iJiii{k) = 7j-^[B+H^^\k7^) + B.Hi'\k7^)] , r/>r/2 (3.4) 

where 

iy=^[a{d-l--f)-np + l] 
^ = ^{niPi+-fi - 1) 

^ = ^(m/32 + 272 - 3) (3.5) 

and C is a normalization constant. The Bogoliubov coefficients are given by 
c±{k) = B±/C, and can be fixed by the four conditions obtained matching iIj 
and if)' aX T] = —r]i and r] = 772- 

The coefficients determined in this "sudden" approximation lead, however, 
to an ultraviolet divergence of the energy density of the produced particles. The 
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reason is that, for modes of comoving frequency k"^ higher than the height of 
the potential barrier, the sudden approximation is no fonger adequate, and the 
mixing coefficients should be computed by replacing the potential step with a 
smooth transition of V{rj). In this way one finds, indeed, that the mixing of 
the modes with k > \V\^ is exponentially suppressed with respect to the other 
modes^'^°'^^, and the ultraviolet divergence is avoided. In this paper, however, we 
are mainly interested in the general behaviour of the spectral distribution, and not 
in the details of the transition regime. We shall completely neglect, therefore, the 
frequency mixing of modes which never "hit" the potential barrier, by putting, 
for such modes, c+(A;) ~ 1, c_(A;) ~ 0. This replaces the exponential decay of the 
high frequency side of the spectrum with a cutoff, at an appropriate frequency 
k ~ \V\ ^ . 

Our potential barrier (3.3) has two steps, which satisfy V(r]i) ~ -q^^ » r/J^ ~ 
V{r]2) (for realistic values of the parameters). The propagation of modes with 
r]2^ < k < r]^^ will thus be affected, in our approximation, only by the first back- 
ground transition at 77 = ?7i. In this frequency band, the Bogoliubov coefficients 
are then defined by c± = A±/C; by matching t/;/, t/;// and their first derivatives at 
rj = rji, and by using the small argument limit of the Hankel functions, we obtain 
(for kr]i < 1) 

1 r Jem-, hri-, 

(3.6) 



1 



2 

(here 7 = r(//)/r(z^), where F is the Euler function, and we have supposed > 0, 
f > when performing the — > limit). 

These coefficients satisfy correctly the Bogoliubov normalization condition, 
|c_|_p — |c_p = 1, and have been obtained in a more particular case-^^, and also 
with a different procedure^'^^, in previous papers. For krji < 1, we shall keep 
the dominant term only, ignoring corrections to the sudden approximation near 
th maximum frequency ki = rj^^, and neglecting also numerical factors of order 
unity, which depend on the model of background evolution (continuity of the scale 
factors and of the dilaton at the transition time), and which do not affect the 
qualitative behaviour of the spectrum. In the rest of the paper, therefore, we shall 
use the expression 

|c_(A;)| = (/cr/i)-!'^-'^! , k^ < k < ki (3.7) 

where A;2 = l/?72 is the frequency corresponding to the height of the barrier V{r]2)- 
Lower frequency modes, /c < /c2, are affected also by the second background 
transition, ai rj = r]2, from the radiation to the matter-dominated regime^'^'^. In 
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this frequency sector the Bogohubov coefficients are given by c± = B±/C, and 
the matching condition provide, for k'r]2 < 1, 



(3.8) 



where 72 = r(cr)/r(/i) for > and a > 0. It may be useful to note that the 
expression can be easily generalized , by performing the product of n Bogoliubov 
transformations, to the case of n background transitions, at i] = r]i, between the 
mode solutions Hi^. and H,y._^-^, with i = 1,2, ...,n. One finds, in generaP^, 



i=l 



7.(^)----±7r^(^)^-^- 



(3.9) 



where 7^ are numerical factors of order unity, and N* = is an overall constant 
phase factor. 

In order to keep only the dominant term of eq.(3.8), for /c < /c2 <^ /ci, we 
have to note first of all that the phenomenological constraints on the time varia- 
tion of the fundamental constant (including G), during the matter and radiation- 
dominated era, imply cr — /i < (see Sec. 5). If ^ — u < (as seems to be indeed the 
case for all the appropriate models of background evolution, see Sec. 6), the second 
term on the r.h.s. of cq.(3.8) is the dominant one. If, on the contrary, fi — 1/ > 0, 
then the first term is dominant (for realistic values of r]i and 772). We shall thus 
use, for the graviton production at low frequencies, 

|c_(/c)| ~ (/cr7i)-l'^-^l(A;r72)^'''"''' , ko < k < k2 (3.10) 

where the — (+) sign refers to /U — z/ < (> 0), and ko is the minimal amplified 
frequency^'^ emerging to-day from the barrier (otherwise stated: crossing to-day 
the Hubble radius Hq^), namely ko = ao-ffo- 

The final number of produced gravitons, for each mode k, is given by |c_(A;)p. 
The corresponding energy density Pg, in the proper frequency interval du, is ob- 
tained by summing over the two polarization states, and is related to c_ by^'^ 

du 



dpg = 2ij\c-\H7ru;^j—^ (3.11) 



The spectral energy density p(a;) = codpg/dw, which is the variable usually 
adopted^'^"^ to characterize the graviton energy distribution, turns out then to 
be parametrized as follows 

p{u;) ~ a;^(A;r7i)-2|'^-^l , k2 < k < ki 
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p(cc;) ~a;^(A;r/i)-2|'^-'^l(/er72)^^'^"''' , ko < k < k2 (3.12) 

For later comparison with present observational data , it is convenient to 
replace all comoving frequencies k by the associated proper frequency u> = k/a{t), 
and to express the spectral distribution in terms of the final curvature scale Hi = 
H^rji), reached at the end of the inflationary phase, Hi ~ (ai?7i)~^ = loi. Since in 
our model rji is also the beginning of the radiation-dominated evolution for a(t), 
it follows that H can be expressed in terms of the radiation energy density p^, as 

H?^[tiy^G,,M , G,^(t) ^ (3.13) 

Note that we have used the Newton constant G ~ M~'^ as the effective gravi- 
tational coupling during the post-inflationary cosmological evolution; the allowed 
deviations from this value turn out to be indeed negligible for our determination 
of the spectral behaviour (see Sec. 5). 

By using eq.(3.13), and by measuring p{u>) in units of critical energy density 
Pc, the spectral distribution (3.12) can be recast finally in the convenient form 
{Q{uj) = p{uj)/pc) 

Q{u!,t) GHfQ^{t) I j , U!2<U!<U!l 

n{u;,t)c^GHfn^{t) ( — ] l—j , uJo<uj<u;2 (3.14) 

where Qj{t) = pj{t)/pc is the fraction of critical energy density present in the 
form of radiation, at the given observation time t. 

This spectrum is parametrized by the scale Hi, and by the kinematical indices 
/i, z/, o", which determine its frequency behaviour. It may be interesting to note that 
the high frequency part of the spectrum is decreasing, flat or increasing depending 
on whether \fi — u\ is larger, equal or smaller than 2. For a primordial phase 
corresponding to isotropic inflation of d = 3 spatial dimensions, with frozen dilaton 
and internal radius (/3 = /3i = 7 = 71 = 0), eq.(3.4) gives, in particular, |/U — z/| = 
1 + a, so that the behaviour of the spectrum is the same as that of the curvature 
scale. For a de Sitter phase (a = 1) one recovers indeed the well known flat 
spectrum^'^^ (O ~ const.), while for superinflation (0 < a < 1) one obtain the 
growing spectrum recently discussed in Ref.22. 

In the general case in which (i 7^ 3, and the additional contributions of a 
dilaton variation (as well as those of dimensional reduction) are included, how- 
ever, the spectral behaviour may be flat or decreasing even if the curvature is 
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growing. What is important to stress is that, in any case, all observational data 
and constraints on the present background of cosmic gravitational waves can be 
translated, thanks to eq.(3.13), into direct information on the curvature scale Hi 
(marking the transition from the primordial inflationary phase to the standard 
decelerated scenario), and on the kinematics of the background evolution. This 
possibility will be discussed in Sec. 5. 

We conclude this section with an estimate of the transition frequencies cji and 
UJ2. At the present time to, the minimal proper frequency ujq is determined by the 
to-day value of the Hubble radius, i.e. ujq = Hq 10~^^Hertz. The frequency 
u!2, corresponding to the matter-radiation transition, can be easily related to uq 
by noting that a{t) ~ fs during the matter-dominated regime, so that 



u^^k2^ H2a2 ^ (t_oy ^ faoy 
u)o ko Hotto \t2j V«2y 

On the other hand, the radiation temperature evolves adiabatically (aT = const.), 
so that the ratio (3.15) can be expressed in terms of the temperature T2 at the 
transition time, 

1 

1^2 { T2 



- W (3.16) 

"-1^0 \ / 

where Tq ~ l^K is the present temperature of the radiation background. 

In a similar way we can relate ujq to the maximal cutoff frequency u)i, which 
depends on the final curvature scale Hi. We can put, in fact, 

LOi ki Hiai f Hiai\ f H2a2\ (3 17) 



u}Q ko HottQ \H2a2 J \ Hotto J 

and we note that, during the radiation dominated evolution, a ~ ~ H~^. 
We have, moreover, H2 ~ lO^i^o and (in units of Planck mass) Ho ~ 10~^^Mp; 
therefore 

^1 . -,^2fHiY fMpY fHoY ,^29 fH 



4. The squeezing parameter 

Another phenomenological signature of the primordial cosmological transi- 
tions, encoded into the cosmic gravity-wave background, is the squeezing param- 
eter which characterizes the quantum state of the gravitons produced from the 
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vacuum^^. This parameter is directly related to the Bogoliubov coefficients, and 
is thus sensible to all the various components of the production process, includ- 
ing a possible variation of the dilaton background, just like the spectral energy- 
distribution. 

The graviton production discussed in the previous section is based on the 
expansion of the gravitational perturbation in terms of \in > and \out > states, 
namely 

ip{k,r]) = bipin + b^*^ 

%l){k, T]) = a%l)out + oVout (4-1) 

for each mode k. The two sets of solutions are connected by a Bogoliubov trans- 
formation which, when expressed in terms of the "in" and "out" mode solutions, 
takes the form 

/ in- \ In. n \ f on . \ 

(4.2) 



'>Pin \ _ ( Cj^ C- \ ( ''Pout 



where c± are defined, according to eq,(3.4), as c± = B±/C. The equivalent rela- 
tion among the corresponding annihilation and creation operators of the second- 
quantization formalism is then 

If the Bogoliubov transformation is parametrized by two real numbers, r and 
$, in such a way that 

c+ = cosh r , c*_ = -e^*^ sinh r (4.4) 

the transformation (4.3) can be rewritten as 

a = S^bS , a"l" = S^b^S (4.5) 

where (5 is a unitary operator defined by 

5 = exp[^z(6t)2_ , z = re'''^ (4.6) 

This is a so-called "squeezing" operator: when applied to the vacuum (or, more 
generally, to a coherent state), generates a state for which the quantum fiuctuations 
of the operator X ~ 6 + 6^ (or its canonical conjugate) can be arbitrarily squeezed 
for a suitable choice of r (see for instance Ref.(24)). In particular, AX —>■ for 
r —>■ oo. 
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The cosmic gravitons arising from the background transitions are thus pro- 
duced in a squeezed state, with a parameter r which, according to eq (4.4), is given 

by 

r = ln(|c_| + V|c_|2 + 1) (4.7) 

According to the model of background evolution considered in the previous section, 
and for uj > uj2, the relic graviton background may be characterised, in general, 
by the following squeezing parameter 

r{oj) ~ In |c_| ~ — l/i — i/| In [ — 



— Ia* ~ ^1 



"^1 / / s 

(4.8) 



25-lnf^) + llnr^'| 
\HertzJ 2 \MpJ 



(we have used eq.(3.7) for c_, and the estimate (3.18) for uji). 

The first term in eq.(4.8) is expected to be the dominant one, at least in 
the range of frequencies accessible, in a (hopefully) not too distant future, to a 
direct observation^'^. The second term takes into account the variation of r with 
frequency, and the third term provides a correction if the transition curvature 
scale is different from the Planck scale. A direct measurement of this parameter, at 
some definite value of frequency, would provide then significative information both 
on the curvature scale Hi, and on the background (dilaton included) dynamical 
evolution, through the |a* ~ ^^1 dependence. 



5.PhenomenoIogical constraints on the graviton spectrum 

The present energy distribution of a cosmic gravity- wave background is mainly 
constrained by three kind of direct observations"^'"^: the absence of fluctuations in 
the millisecond pulsar-timing data, the critical density value, and the isotropy of 
the cosmic microwave background radiation (CMBR). The flrst one applies on a 
narrow frequency interval around ~ 10~^Hertz, while the other two at all 
frequencies (the third one provides a bound which is frequency-dependent). Their 
relative importance, and the frequency at which they provide the must significant 
constraint, depend on the slope of the graviton energy spectrum Q{uj). 

For a stochastic graviton background, the bound on the spectrum following 
from the CMBR isotropy constrains the wave amplitude h(uj), and scales like u~'^. 
It provides then the most significative bound at the minimum frequency ujq (where 
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it implies^ f2 < 10~^), unless we have a spectrum which in its low frequency band 
(a; < (JJ2) grows faster than a;^. If the spectrum is growing at all frequencies, 
however, the most significant constraint, for the present values of experimental 
data, is provided in any case by the critical density bound ~ 1, applied to the 
highest frequency ui . 

According to our three component model of background evolution, the spec- 
trum may be increasing at low frequencies, and simultaneously flat or decreasing 
in the high frequency sector, only if (see eq.(3.14)) 

II , 2 < I// - < 2 + |cr - //| (5.1) 

Even in such a particular case, however, the growth of the low frequency sector 
cannot be significantly faster than cu^, since, as we shall see later, \a — ii\ is not 
allowed to be notably larger than 1 by the present limits on the variation of 
the fundamentals constants. Therefore, the energy distribution of the graviton 
background can be significantly constrained by imposing on eq.(3.14) the three 
following bounds: 

Q{oji) < Qc , ^{i^p) < , ^{(^0) < (5.2) 

where Up ~ 10~^Hertz, and Qc, ^p, are the present value of the bounds on the 
energy density imposed, respectively, by critical density, pulsar timing data, and 
CMBR isotropy. 

For our discussion of the constraints, it may be convenient to simplify the 
notations by defining the variables 

x=\n-u\ , y=\a-ii\ , ^ = log(^^) (5-3) 

By using fl^{to) ~ 10~^ for the present critical fraction of radiation energy density, 
and by inserting in eq.(5.2) the values of ujq, uji, uj2 determined in Sec. 3, the 
three constraint equations in the parameter space, for our model of background 
evolution, can thus be written 

z <2+ ^log^c 
z < -(SO + logfip) -38 

X 

2 < - (120 T % + log Qi) - 58 (5.4) 

X 
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They follow, respectively, from the critical density, pulsars and isotropy bounds, 
and they define an allowed region in the {x,y,z) space which provides information 
on the past evolution of our Universe. 

In order to discuss the extension of this region it should be noted, first of all, 
that the range of variation of the variable y, 

y=\a-fi\ = ^|ni(/?2 - A) + 272 - 71 - 2| (5.5) 

which parametrizes the time evolution of the dilaton and of the compactification 
radius during the matter and radiation dominated era (recall eq.(3.2)), is severely 
constrained by the present bounds on the variation of the fundamental constants. 

Indeed, in a Brans-Dicke frame, and in a higher dimensional context with 
n — D — 4 dimensions lying in a compact internal space, with scale factor b{t), the 
effective four dimensional Newton constant Gn evolves in time like Gjv ~ /b^ . 
We have then 

^ = i-n\ (5.6) 

Gat 

During the matter dominated era the variation of the extra spatial dimensions is 
constrained by^^ 

\h/h\ < IQ-^Ho (5.7) 

and the variation of Gat by^^ 

\Gn/Gn\ < (5.8) 

where we have taken for Hq the largest value allowed to-day, Hq ~ 10~^^yr~^ . 
These two bounds imply \(p\ < 10~^Hq. But, according to our parametrization 
(3.2), ^ = 72 and b/b = -P2H/2. It follows that 

I/32I < 10"^ , I72I < 10"' (5.9) 

Consider now the radiation-dominated era. During this phase, the best lim- 
its on <j) and b are obtained from the primordial nucleosynthesis. Denoting by 
bnuch Gnuci-i ct^d by 60 : Gq, the values of the radius of the internal space and of 
the Newton constant, at the epoch of nucleosynthesis and at the present epoch, 
respectively, one obtains that the change of b must be bounded by^^'^^ 

%^ = l + e , |e|<10-2 (5.10) 
bo 
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while the change of G is constrained by^^ 

= l + e , |e|<3xl0-^ (5.11) 



Crnucl 

Go 



Translated into limits on the time variation of b and (f), according to the 
parametrization (3.2), they imply 

|/3i|<10-3 , |7i|<10-i (5.12) 

The dilaton contribution is thus the dominant source of uncertainty in the 
value of the parameter y. Even taking into account the maximum allowed uncer- 
tainty, however, it follows from eqs.(5.9) and (5.12) that 

0.9<y<l.l (5.13) 

A first rough evaluation of the allowed region in the parameter space is thus 
obtained by fixing y = 1 in eqs.(5.4) (the allowed deviation of y from 1 is too 
small to be significant in view of our previous approximations). 

We have to insert, moreover, in eq.(5.4) the values of the bounds implied by 
the present experimental data. We shall put Oc = 1 (in order to avoid that the 
produced gravitons over-close our present universe), Qp = 10~^ as implied (at 
the 99% confidence level) by recent results from pulsar timing^^, and Qi = 10~^, 
following from the constraint h < 10~^ on the gravity-wave amplitude. With 
this data, the constraint eqs.(5.4) become 

z<2 

z< 38 

X 

;z < i(112T4) -58 (5.14) 

We recall that the negative (positive) sign in the last equation corresponds to 
H < u {p > u) . It should be mentioned, moreover, that in the context of a more 
stringent analysis, the first bound z < 2 could be replaced by z < 0, following from 
the fact that early nucleosynthesis seems to imply^^ , at high frequency, O < 10~^ 
for the energy density distribution of massless particles. This would correspond to 
a maximum scale Hi < Mp instead of lO^Mp. This conclusion is, however, model 
dependent, and in this paper we prefer to rely on constraints following directly 
from observations. 
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The allowed region of the {x,z) plane delimited by eqs.(5.14) is illustrated in 
Fig.l. Because of the uncertainty of the experimental data, which has not been 
completely taken into account in our discussion, and because of the approximations 
made, this figure is expected to give only a qualitative picture of the phenomeno- 
logical scenario. Nevertheless, we can draw from our analysis the following general 
conclusions. 

1) There is a maximum allowed value for the curvature scale Hi at the epoch 
of the transition from the phase of accelerated expansion, dilaton growth and 
dimensional reduction, to the decelerated radiation-driven evolution, i.e. Hi < 
lO^Mp. 

2) Models characterised by a sufficiently high scale. Hi > 10~^Mp , are con- 
strained by pulsar timing if fx > v, and by CMBR isotropy if fx < i/. 

3) For any given scale Hi lower than the maximum one there is a limiting slope 
of the spectrum, below which that scale is forbidden. Within our approximations, 
the limiting slope for a scale Hi is fixed by 

X < r- (5.15) 

58 + log (i^) 

if H < u, and 

74 

X < ^ (5.16) 

38 + log (i;) 

if fi > u. In the first case (which corresponds to all the physical models considered 
in the next section), the maximum scale lO^Mp is allowed for x ~ 1.8, while the 
Planck scale can be reached for x ~ 54/29 ~ 1.86. A four dimensional inflationary 
background, with frozen dilaton and radius of the internal dimensions (7 = /3 = 0, 
d — 3), corresponds in particular to n — u — —a — 1 < 0, and the Planck scale is 
thus reached for a ~ 25/29, in agreement with the results of a previous analysis^^. 

4) Finally, models corresponding to a spectrum which is flat or decreasing 
at high frequencies (i.e. with x > 2), are characterized by a maximum allowed 
scale Hi < IQ-^Mp. We thus recover the well known bound on the scale of a 
four dimensional de Sitter inflation^'^^, since in that case x — \a + 1\ — 2 and one 
obtains the usual flat spectrum. 
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6. String cosmology pre-big-bang 
and other higher-dimensional models 



In the standard cosmological model the curvature is monotonically increasing 
as we go back in time, and blows up at the initial singularity. A possible classical 
alternative to the singularity would seem to be provided by an initial inflationary de 
Sitter phase, at constant curvature, which extends in time indeflnitely toward the 
past. However, as discussed in a recent paper eternal exponential expansion, 
with no beginning, is impossible in the context of the conventional inflationary 
scenario, so that a primordial phase of constant curvature does not help to solve the 
problem of the initial singularity. Moreover, according to the constraints reported 
in the previous section, the constant value of the curvature during the initial de 
Sitter phase should lie at least four orders of magnitude below the Planck scale; 
this may seem unnatural, if one believes that the growth of the curvature is stopped 
and that the primordial curvature becomes stable just because of quantum effects. 

A different alternative has been recently suggested, on the grounds of string 
theory motivations^'^'^^''^'^, in which the singularity is avoided because the cur- 
vature grows up to a maximum (Planckian) value and then decreases back to 
zero. The standard radiation-dominated phase is then preceeded in time by a 
phase with "dual" dynamical behaviour (the curvature and the dilaton are grow- 
ing, H > 0, (p > 0, the evolution is accelerated, a > 0), called^^ "pre-big-bang". 
Particular examples of such a scenario are thus provided also by earlier models of 
superinflation and dynamical dimensional reduction, discussed in the context of 
Kaluza-Klein cosmology^^"^^. 

In this section we want to stress that if the initial conflguration of our model 
of background evolution (i.e. for i] < — ?7i, see eq(3.11)) corresponds to a pre-big- 
bang scenario of this type, the consequent graviton spectrum is always growing 
fast enough to avoid the de Sitter bound H < 10~^Mp (i.e. x < 2), and to allow 
the Universe to inflate up to the maximal curvature scale, consistently with the 
bounds of the previous section. 

Consider indeed the perfect-fluid dominated model of Refs.(34) and (35), 
which describes superinflation and dimensional decoupling, and belongs to the 
class parametrized by eq.(3.1) with 7 = 0. One flnds, for this model, n ~ u = 
and X = \fi — iy\ = 0.5. The model of Ref.(36) (based on the toroidal compact- 
iflcation of D=ll supergravity), corresponds to 7 = 0, a = 0.26 and (3 = 0.22, 
and gives ji — v — —0.49 < 0. The model of string-driven inflation of Ref.(33) 
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has 7 = 0, n > 10, and for d = 3 it gives ji — v = (4 — n)/3n < 0. Finally, a 
typical pre-big-bang model^^, dual to the standard radiation phase, satisfies the 
Brans-Dicke equations (2.2) and (2.3) (with a; = 1) for 



and implies 



3 + a + n 



II- u = <0 (6.2) 

^ 3+d+n ^ ' 



For all these models we have fi — u < 0, and |/U — z/| < 1.8 (for any allowed number 
of internal dimensions), so that their final curvature scale is only constrained by 
the closure density bound. 

We want to comment, finally, on the possibility that the CMBR anisotropy 
recently measured"^^ by COBE be at least partially determined, at the quadrupole 
level, by a cosmic graviton background. It has been already pointed out'^^, indeed, 
that a stochastic background of gravitational waves with flat spectrum, generated 
by a primordial de Sitter inflationary phase, could produce the entire observed 
signal, provided de Sitter inflation occurred at a vacuum energy scale Mpv^^^ ~ 
1.5 X 10^^ GeV (at the 95% confldence level). This translates into a value of the 
Hubble constant H = {8nMp/3)v^/'^ ~ 10~^Mp, which is not in conflict with the 
previously reported bound {Hi ~ 10""^ for x = 2, see Fig.l). 

It should be noted, however, that a four-dimensional de Sitter inflation is not 
the only primordial phase which can be associated to a flat graviton spectrum. 
Indeed, in a more general higher-dimensional Brans-Dicke scenario, all the models 
with l/i — z/| = 2 provide a flat high frequency spectrum. Included in this class, in 
particular, are all the {d+ 1) -dimensional models providing a phase with variable 
dilaton and isotropic superinflationary expansion, characterized in conformal time 
(according to eq.(3.1)) by the power a = 2/ {d — 1 — 'j). 

It remains still open, moreover, the interesting possibility that the COBE 
anisotropy may be fltted a by non-flat graviton spectrum^^ with, in particular, 
x < 2, as predicted by the string pre-big-bang models. In this case we may 
expect, according to Fig.l, that the COBE data will select a higher transition 
scale Hi, and in such a case they could be interpreted, instead of a flrst direct 
evidence, via gravitational waves, for the GUT scenario^^, as evidence for the 
dilaton-driven string cosmology scenario. In order to discriminate between these 
two (exciting) possible interpretations, however, one should try to probe directly 
the energy density of the cosmic graviton background at some given frequency. 
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for example through a gravity wave detector (such as LIGO^), or by means of 
astrophysical methods (such as timing measurements of miUisecond pulsars^^). 



7. Conclusions 

In this paper we have considered a three-component model of cosmological 
evolution in which the standard radiation and matter-dominated expansion of the 
three-dimensional space is preceeded in time by a general d-dimensional phase of 
accelerated (i.e. inflationary) expansion. We have included, moreover, a possible 
variation of the effective gravitational coupling and of the compactification scale, 
parametrized, respectively by a logarithmic time dependence of the dilaton field, 
and by a power law evolution of the internal scale factor. 

We have shown that the linearised equation for a metric fluctuation, obtained 
by perturbing the Brans-Dicke equations around this background, contains a cou- 
pling of the perturbation to the background metric and to the dilaton field 0. 
As a consequence, both the dimensional reduction process and the variation of 
G (via (j)) contribute (besides inflation) to the process of the amplification of the 
gravitational perturbations (i.e. to the graviton production). 

We have computed the spectral distribution Q,{uj) of the energy density stored 
to-day in a cosmic graviton background (and the associated squeezing parameter 
r(a;)), taking into account all possible contributions. The frequency behaviour 
of the spectrum turns out to be clearly related to the temporal behaviour of the 
background fields {gf^y and 0); the observational constraints on Vt{u3) provide then 
significative information both on the kinematics of the background evolution, and 
on the curvature scale Hi characterizing the transition from the primordial in- 
flationary phase (with variable dilaton), and the standard radiation-dominated 
phase. 

We have shown, in particular, that for flat or decreasing spectra the transition 
scale cannot overcome a maximum value which lies, typically, four orders of mag- 
nitude below the Planck scale. For growing spectra, on the contrary, the allowed 
transition scale can be as high as the Planck one (and somewhat higher). 

We have stressed, finally, that the contribution of the dilaton background to 
the cosmic production of gravitons may simulate the usual flat four-dimensional de 
Sitter spectrum, even if the inflationary evolution of the scale factor is not of the 
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exponential type, and the curvature scale is growing, instead of constant, during 
the inflation. As a consequence, one could try to interpret the recently measured 
COBE anisotropy not only as evidence for de Sitter inflation at the GUT scale^^, 
but also (alternatively) as a possible evidence for a dilaton-driven string cosmology 
scenario-'^^'-'^^. 
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Figure caption 



Fig.l 

The maximum allowed value of the transition scale Hi (in units of Planck mass), 
versus the parameters determining the kinematics of the background evolution. 
The allowed region with Hi ~ 10^ Mp extends from x — 1.8 down to x = 0. 
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